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This  report  describes  the  variability  of  the  microwave  (3  to  30  GHz)  circuit 
design  parameter,  effective  earth  radius  factor,  at  two  sites  which  served  as 
test  cases:  Balboa,  Panama,  and  Wiesbaden,  Germany.  Median  effective  earth 
radius  factors,  K,  derived  from  meteorological  data  for  the  first  100  m altitude 
at  these  two  sites,  were  1.32  for  Wiesbaden  and  1.68  for  Balboa.  Median  K fac- 
tors for  84  other  world  localities  are  included  in  table  1.  Numerous  medians 
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and  Wiesbaden  is  illustrated  with  sets  of  equivalent  earth  profile  curves  en- 
closing 68.3,  95.4,  99.7,  and  100  percent  of  the  data  and  by  K factor  cumulative 
distribution  functions  bounded  by  respective  90  and  99  percent  confidence  limits 
Effective  earth  radius  factors  were  highly  variable,  particularly  at  Balboa, 
where  the  meteorological  variability  is  high  within  the  first  100  m altitude. 
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INTRODUCTION 


The  US  Army  Communications  Command  (USACC),  under  the  provisions  of  AR 
10-13,  is  responsible  for  performing  electromagnetic  wave  propagation 
engineering  services  for  agencies  of  the  US  Army  and  other  government 
offices  and  for  maintaining  the  operational  electromagnetic  compatibility 
(EMC)  program  area  [lj.  These  responsibilities  have  been  functionally 
assigned  to  the  Electromagnetic  Engineering  Office  (EMEO)  located  within 
the  US  Army  Communications-Electronics  Engineering  Installation  Agency. 
Under  this  mandate  one  of  the  tasks  the  EMEO  performs  is  to  provide  en- 
gineering consultation  and  RF  system  performance  analysis  and  design  for 
line-of-sight  (LOS)  super  high  frequency  (SHF,  3 to  30  GHz)  communications 
links  [1],  For  each  SHF  link,  transmitter  and  receiver  spacing,  antenna 
sizes  and  heights,  an  all-year  median  path  attenuation  figure,  and  a per- 
centage of  time  during  which  a specified  service  should  be  achieved  are 
all  determined  as  part  of  the  circuit  design.  To  determine  the  above 
quantities,  an  average  radio  wave  refractivity  gradient  is  usually  as- 
sumed, converted  to  an  effective  earth  curvature  correction,  and  entered 
into  the  calculations.  Tropospheric  radio  wave  refraction  gradients, 
however,  are  a function  of  atmospheric  relative  humidity,  temperature, 
and  pressure,  three  highly  variable  parameters.  Recognizing  the  vari- 
ability of  refractivity  gradients  as  a function  of  time  and  world  local- 
ity, the  EMEO  requested  that  the  Atmospheric  Sciences  Laboratory  (ASL) 
perform  the  statistics  which  delineate  the  range  of  effective  earth  ra- 
dius values  to  be  encountered  at  specified  world  localities  by  electro- 
magnetic waves  used  in  SHF  systems  [2]. 

This  report  describes  the  mathematical  procedure  used  in  interpreting 
the  meteorological  data  and  the  statistical  technique  employed  to  derive 
effective  earth  curvature  variability  and  presents  the  results  obtained 
for  Balboa,  Panama,  and  Wiesbaden,  Germany.  The  fundamental  ideas  of 
Samson  [3]  in  handling  refractivity  gradients  were  studied;  then  his 
original  Wiesbaden  and  Balboa  meteorological  data  were  obtained  and  pro- 
cessed. Processing  of  the  meteorological  data  involved  the  determina- 
tion of  over  4700  refractivity  gradients  at  each  station.  Sincp  most 
LOS  microwave  transmitters,  repeaters,  and  receivers  employ  high-gain 
directional  antennas  located  within  the  first  100  m above  ground  level, 
all  atmospheric  refractivity  gradients  derived  here  are  for  that  inter- 
val between  the  earth's  surface  and  100  m altitude.  Effective  earth  ra- 
dius K values  were  determined  corresponding  to  the  median  value  of  the 
refractivity  gradient  at  each  station  and  to  selected  percentile  devia- 
tions of  the  data  from  the  median.  To  provide  US  Army  microwave  engi- 
neers with  a design  tool,  equivalent  earth  profiles  were  plotted  for 
these  K values.  Cumulative  distribution  functions  of  K and  the  90  and 
99  percent  confidence  limits  for  each  K value  were  also  determined  from 
the  data.  Details  of  the  aforementioned  analysis  are  included  as  part 
of  this  report. 
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METEOROLOGICAL  AND  REFRACTIVITY  DATA 


The  primary  sources  of  meteorological  data,  excluding  surface  obser- 
vations, are  those  measurements  of  relative  humidity,  pressure,  and 
temperature  obtained  with  standard  balloon-borne  radiosondes  in  use 
throughout  the  world.  Radiosonde  observations  of  temperature,  pressure, 
and  relative  humidity  are  reported  at  fixed  atmospheric  pressure  levels 
termed  "mandatory  levels"  and  at  "significant  levels."  Significant 
levels  are  those  where  an  appreciable  change  in  temperature  or  relative 
humidity  as  a function  of  altitude  occurred  before  the  balloon  reached 
the  next  mandatory  level.  Consequently,  mandatory  levels  differ  in 
altitude  because  of  differences  in  terrain  elevation  at  each  sounding 
station.  Significant  levels  differ  in  altitude  due  to  the  variable  tem- 
perature lapse  rates  encountered  during  the  year.  Normally  the  first 
mandatory  level  or  significant  level  above  the  surface  is  at  an  alti- 
tude greater  than  100  m.  Therefore  to  calculate  refractivity  at  100 
m altitude,  and  subsequently  the  refractivity  gradient  for  that  interval, 
an  interpolation  scheme  had  to  be  employed  which  utilizes  the  calculated 
refractivity  values  obtained  from  the  meteorological  measurements  at 
the  surface  and  at  the  lowest  mandatory  level. 

The  expression  for  radio  wave  refractivity  in  the  lower  troposphere,  for 
frequencies  up  to  30  GHz,  is  well  established  [4]  and  is  of  the  form: 

N = 77.6  y + 3.73  x 105  1 (1) 

1 T 

where  P is  the  total  atmospheric  pressure  in  miWibars,  T is  absolute 
temperature  in  degrees  Kelvin,  and  e is  the  partial  pressure  of  water 
vapor  in  millibars.  N is  in  units  of  refractivity,  related  to  the 
usual  index  of  refraction  n by  N = 106(n  - 1).  In  equation  1,  the 
quantity  77.6  P/T  is  normally  referred  to  as  the  "dry  term,"  while 
the  quantity  3.73  x 105e/T2  is  called  the  "wet  term."  Accordingly, 
the  refractivity  may  be  expressed  as 

N = Nd  + Nw  , (2) 


where 

Nd  = 77.6  P/T  , (3) 

Nw  = 3.73  x 105e/T2  . (4) 

Variations  in  refractivity  as  a function  of  altitude  have  been  shown  to 
follow  an  exponential  function  [5]  of  the  form 
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Since  the  assembled  meteorological  data  allow  one  to  calculate  ND  and 
N at  the  earth's  surface  and  at  the  first  mandatory  altitude  of  h, , 

W I 

solutions  for  the  parameters  Ap,  Aw,  Bp,  and  Bw  were  obtained  and  sub- 
stituted into  equations  5 and  6.  Subsequently  equations  5 and  6 were 
summed  to  yield  the  following  refractivity  interpolation  equation: 
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W" 

ND*hI 1 


+ W 


’wr*, 


where 

N(h)  = Total  refractivity  at  altitude  h. 


Np(h0)  = Dry  term  refractivity  at  earth's  surface,  hQ, 


Nn(h, ) = Dry  term  refractivity  at  altitude  h,  for  which  data  is 
available, 

N,(h  ) = Wet  term  refractivity  at  earth's  surface,  h . 
w o o 

N (h, ) = Wet  term  refractivity  at  altitude  h,  for  which  data  is 
available. 

Each  interpolation  yielded  a refractivity  magnitude  at  an  altitude  of 
100  m,  which  when  subtracted  from  the  refractivity  magnitude  obtained 
at  the  earth's  surface,  provided  the  desired  refractivity  gradient, 
AN/Ah.  The  above  procedure  shows  that  AN/ Ah  is  assumed  to  be  linear 
in  the  first  100  m of  altitude.  For  the  locality  of  Balboa,  Panama, 

4917  refractivity  gradients  were  computed  with  radiosonde  measurements 
of  the  period  from  1951  to  1956.  During  the  sunrise  atmospheric  heating 
and  sunset  cooling  time  periods,  refractivity  gradients  exhibit  maximum 


variability  and  definite  changes  in  absolute  magnitude.  Thus  in  an 
analysis  of  this  type,  it  is  desirable  to  have  a meteorological  sample 
of  the  sunrise  and  sunset  periods  included.  However,  of  the  total  data 
set,  only  380  samples  were  representative  of  near  sunset  or  sunrise 
hours;  the  remainder  of  the  data  set  consists  of  a balanced  mixture  of 
totally  night  or  daytime  meteorological  data.  For  Wiesbaden,  Germany, 
4798  refractivity  gradients  were  derived  from  data  recorded  during  the 
period  1951  to  1957.  All  Wiesbaden  meteorological  data  are  from  either 
daytime  or  night  radiosonde  measurements;  no  sunrise  or  sunset  samples 
were  available.  Consequently,  the  variability  results  to  be  presented 
here  are  inherently  conservative. 

Ordinal  listings  of  the  100-m  refractivity  gradients  for  Balboa  and 
Wiesbaden  were  compiled  and  the  median  refractivities  of  -64  units/km 
and  -38  units/km,  respectively,  were  obtained.  Initial  attempts  to 
work  with  average  refractivity  gradients  and  deviations  from  the  average 
value  proved  meaningless  because  the  frequency  distribution  of  the  re- 
fractivity gradients  was  found  to  be  non-Gaussian.  Figures  1 and  2 are 
refractivity  gradient  histograms  which  approximate  the  frequency  distri- 
butions at  Balboa  and  Wiesbaden.  These  figures  show  that  the  distribu- 
tion is  non-Gaussian  or  at  least  that  the  sample  was  insufficient  to  a 
Gaussian  distribution.  Therefore,  nonparametric  statistics  were  employed 
to  determine  a median  value  and  percentile  distributions  about  the  me- 
dian. Annotated  on  the  histograms  are  the  median  values  and  percentile 
deviations  from  the  median  obtained  directly  from  the  ordinal  listing. 

The  34th  and  47.5th  percentile  deviations  denote  how  68  and  95  percent 
of  the  data  are  distributed  about  the  median.  For  Balboa,  95  percent 
of  the  gradient  data  is  within  -202  and  -1  units/km;  while  for  Wiesbaden, 
95  percent  of  the  gradients  have  values  ranging  from  -86  to  +16  units/km. 


EFFECTIVE  EARTH  CURVATURES 

Radio  transmission  path  profiles  used  to  design  microwave  communications 
links  must  account  for  the  earth's  curvature  and  radio  beam  refractivity. 
To  expedite  circuit  design  work,  the  radio  beam  refraction  is  normally 
combined  with  the  earth's  curvature  to  define  an  effective  earth  radius 
usually  designated  as  Kr,  V/here  r is  the  true  earth  radius  and  K is  the 
atmospheric  refraction  dependent  modification  factor.  Depiction  of  the 
earth  with  an  effective  earth  radius  of  Kr  allows  the  radio  beam  to  be 
drawn  as  a straight  line  as  shown  in  figures  3a  and  3b.  To  further 
facilitate  path  profile  plotting,  both  curves  can  be  transformed  to  a 
flat  earth  reference  frame  and  a radio  wave  beam  with  curvature  Kr  as 
shown  in  figure  3c  or  3d.  Variations  in  the  K factor,  induced  by  mete- 
orological changes,  may  then  be  depicted  by  several  plots  represent- 
ing the  range  of  probable  Kr  magnitudes. 

In  this  section  the  equations  used  to  derive  all  the  K factors,  the 
median  K,  percentile  deviations  from  the  median  K,  and  equivalent  earth 
profile  curves  as  shown  in  figure  3d  are  described.  Actual  values  and 
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Figure  1.  Balboa,  Panama,  refractivity  gradient  histogram.  Figure  2.  Wiesbaden,  Germany,  re fr activity  gradient  histogram. 


Median  AN/Ah, 


Pth  percentile  value  of  AN/Ah  obtained  oirectly  from  ordinal 
listing 


Equivalent  earth  profile  curves,  like  those  shown  in  figures  3c  and  3d, 
were  constructed  by  using  the  simple  geometric  relationship  [6] 
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where 

h = height  above  path  center  at  horizontal  distance  d, 
d = horizontal  distance  from  path  center, 

K = equivalent  earth  radius  factor. 


Graphs  employing  relationship  shown  in  equation  11  display  the  variabil- 
ity of  radio  beam  refractivity  for  expected  magnitudes  of  K determined 
from  statistical  analysis  using  pertinent  meteorological  data. 

Employing  equation  9,  median  K values  of  1.68  and  1.32  were  calculated 
for  Balboa  and  Wiesbaden,  respectively.  Note  here  that  the  Wiesbaden 
K value  of  1.32  is  in  agreement  with  the  K = 4/3  value  conventionally 
used  throughout  the  world;  however,  the  Balboa  median  differs  signifi- 
cantly from  4/3.  Reference  5 provides  a catalog  of  median  refractivity 
gradients  for  86  world  locations.  These  gradients  have  also  been  con- 
verted to  median  K values  employing  equation  9 and  are  listed  in  alpha- 
betical order  in  table  1. 

Percentile  deviations  from  the  median  equivalent  to  standard  deviation 
percentage  ranges  commonly  used  in  Gaussian  distributions  (a,  2a,  etc.) 
were  identified  in  the  refractivity  gradient  ordinal  listing  and  used  in 
equation  10  to  compute  the  corresponding  K values.  The  median  K values 
and  the  K's  obtained  for  the  68.3th,  95.4th,  99.7th,  and  100th  percentil 


TABLE  1 . EQUIVALENT  EARTH  RADIUS  FACTOR  MEDIANS 


City  State/Country 

Approximate  Geographic 
Coordinates 

Median  K 

Abidjan,  Ivory  Coast 

5 N 

4 W 

K 

s 

1.35 

Aden,  Yemen 

13  N 

45  E 

K 

s 

1.55 

Anchorage,  Alaska,  US 

61  N 

150  W 

K 

* 

1.39 

Aloulef,  Algeria 

27  N 

1 E 

K 

= 

1.35 

Argentia,  Newfoundland 

47  N 

54  W 

K 

= 

1 .32 

Athens,  Greece 

38  N 

24  E 

K 

* 

1.38 

Atlanta,  Georgia,  US 

34  N 

84  W 

K 

= 

1.78 

Balboa,  Canal  Zone 

9 N 

80  W 

K 

= 

1.68 

Bangui,  Central  African 

4 N 

19  E 

K 

= 

1.73 

Republ ic 

Barrow,  Alaska,  US 

71  N 

157  W 

K 

s 

1.33 

Bitburg,  Germany 

50  N 

7 E 

K 

= 

1.34 

Bordeaux,  France 

45  N 

1 E 

K 

= 

1 .62 

Brownsville,  Texas,  US 

26  N 

97  W 

K 

= 

1.65 

Brussel,  Belgium 

51  N 

4 E 

K 

= 

1.54 

Calcutta,  India 

23  N 

88  E 

K 

s 

2.38 

Charleston,  South  Carolina,  US 

33  N 

80  W 

K 

= 

1 .51 

Chitta,  USSR 

52  N 

113  E 

K 

1.37 

Clark  Field,  Luzon, 

15  N 

121  E 

K 

= 

1 .52 

Philippines 

Cocoa  Beach,  Florida,  US 

28  N 

81  W 

K 

= 

1.6 

Columbia,  Missouri,  US 

39  N 

92  W 

K 

= 

1.5 

Coral  Harbor,  Northwest 

64  N 

83  W 

K 

1 .43 

Territories,  Canada 

Curacao,  Netherlands 

12  N 

69  W 

K 

= 

1.67 

West  Indies 

Dakar,  Senegal 

15  N 

17  U 

K 

= 

157.0 

Denver,  Colorado,  US 

40  N 

105  W 

K 

= 

1.34 

El  Paso,  Texas,  US 

32  N 

106  W 

K 

= 

1.44 

Ft.  Lamy,  Chad 

12  N 

15  E 

K 

= 

-78.5 

Ft.  Smith,  Northwest 

61  N 

112  W 

K 

= 

1.35 

Territories,  Canada 

Ft.  Trinquet,  Mauritania 

25  N 

12  W 

K 

= 

1.45 

Gibraltar  (British  Colony) 

36  N 

5 W 

K 

= 

1.6 

Spain 

Gross  Rohrheim,  Germany 

50  N 

8 E 

K 

= 

1.54 

Guantanamo,  Cuba 

20  N 

75  W 

K 

= 

1.45 

Gur'Yey,  Kazah,  USSR 

47  N 

52  E 

K 

= 

1.34 

Hannover,  Germany 

52  N 

10  E 

K 

= 

1.47 

Hilo,  Hawaii,  US 

20  N 

155  W 

K 

= 

1.65 

Hong  Kong  (British 

22  N 

114  E 

K 

= 

1.59 

Colony)  China 

10 

til 
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Table  1 (cont) 


City  State/ Country 

Approximate  Geographic 
Coordinates 

Median  K 

Isachsen,  NM  Territories, 

79 

N 

103 

W 

K 

= 1.45 

Canada 

Joliet,  Illinois,  US 

41 

N 

88 

W 

K 

= 1.51 

Karaganda,  Kazakhstan, 

50 

N 

78 

W 

K 

= 1.33 

USSR 

Key  West,  Florida,  US 

25 

N 

82 

W 

K 

= 1.52 

Khartom,  Sudan 

16 

N 

33 

E 

K 

= 1.34 

Koror  Island,  Pacific 

7 

N 

134 

E 

K 

= 1.94 

La  Coruna,  Spain 

43 

N 

8 

W 

K 

= 1.34 

Long  Beach,  California,  US 

34 

N 

118 

W 

K 

= 1.44 

Long  Xuyen,  Viet  Nam 

10 

N 

105 

E 

K 

= 1.0 

(Only  1 month  of  data 

available;  morning  hours) 

L'Vov,  Ukrain,  USSR 

50 

N 

24 

E 

K 

= 1.4 

Madrid,  Spain 

40 

N 

4 

W 

K 

= 1.35 

Majuro,  Marshall  Islands, 

7 

N 

171 

E 

K 

= 1.64 

Pacific 

Mazatlan,  Mexico 

23 

N 

106 

W 

K 

= 1.8 

Miami , Florida,  US 

26 

N 

80 

W 

K 

= 1.65 

Midway  Island,  Pacific 

28 

N 

177 

W 

K 

= 1.54 

Moscow,  USSR 

56 

N 

38 

E 

K 

= 1.34 

Muharraq,  Bahrain 

26 

N 

51 

E 

K 

= -31  .< 

(Persian  Gulf) 

New  Delhi,  India 

29 

N 

77 

E 

K 

= 2.57 

New  York,  New  York,  US 

41 

N 

74 

W 

K 

= 1.34 

Niamay,  Niger 

13 

N 

2 

E 

K 

= 5.81 

Nicosia,  Cyprus 

35 

N 

33 

E 

K 

= 1.32 

Oakland,  California,  US 

38 

N 

122 

W 

K 

= 1.48 

Odessa,  Ukraine,  USSR 

46 

N 

31 

E 

K 

= 1.4 

Okhotsk,  Khabarovsk 

59 

N 

143 

E 

K 

= 1.38 

Territory,  USSR 

Ostersund,  Sweden 

63 

N 

15 

E 

K 

= 1 .33 

Palma,  Majorca 

40 

N 

3 

E 

K 

= 1.29 

Port  Lyautey,  Morocco 

34 

N 

7 

W 

K 

= 1.6 

Saigon,  Viet  Nam 

11 

N 

107 

E 

K 

= 1.8 

San  Diego,  California,  US 

33 

N 

117 

W 

K 

= 1.5 

San  Juan,  Puerto  Rico 

18 

N 

66 

W 

K 

= 1.8 

Seattle,  Washington,  US 

47 

N 

122 

W 

K 

= 1.41 

Shemya,  Alaska,  US 

53 

N 

174 

E 

K 

= 1.3 

Ship  "C"  (Atlantic  Ocean) 

53 

N 

35 

W 

K 

= 1.34 

Ship  "K"  (Atlantic  Ocean) 

45 

N 

16 

W 

K 

= 1.41 

Ship  "M"  (Norwegian  Sea) 

66 

N 

2 

E 

K 

= 1.34 

Ship  "V"  (Pacific  Ocean) 

34 

N 

164 

E 

K 

= 1.65 

Singapore  (South  of  Malay 

1 

N 

104 

E 

K 

= 1.99 

Peninsula) 

Stuttgart,  Germany 

49 

N 

9 

E 

K 

= 1.47 

Swan  Island,  West  Indies 

17 

N 

84 

W 

K 

= 1.62 

11 


Table  1 (cont) 

City  State/Country 

Approximate  Geographic 
Coordinates 

Median  K 

Syktyvkay,  Komi,  USSR 

62  N 

51  E 

K = 1.33 

Tampa,  Florida,  US 

28  N 

83  W 

K = 1 .65 

Tan  An  Viet  Nam 

11  N 

106  E 

K = 1 .8 

Tashkent,  U Bekistan,  USSR 

41  N 

69  E 

K = 1 .47 

Tateno,  Japan 

36  N 

140  E 

K = 1 .48 

Tatoosh  Island,  Washington,  US 

48  N 

125  W 

K = 1.43 

Tura,  Krasnoyarsk 

64  N 

100  E 

K = 1.45 

Territory,  USSR 

Verkhoyansk,  North  Yakut, 

67  N 

133  E 

K = 1 .59 

USSR 

Vladivostok,  USSR 

43  N 

133  E 

K = 1 .39 

Wake  Island,  Pacific 

19  N 

167  E 

K = 1 .85 

Washington,  DC,  US 

39  N 

77  W 

K = 1.34 

Wiesbaden,  Germany 

50  N 

8 E 

K = 1 .32 

(Federal  Republic) 
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regions  about  the  Balboa  and  Wiesbaden  median,  respectively,  are  tabu- 
lated in  figures  4 and  5.  Inspection  of  the  K column  in  figures  4 and 
5 reveals  three  ranges  of  the  K magnitudes  which  are  of  interest:  0 < 
k < 1,  1 ^ K < °°,  and  K negative. 

These  K data  sets  indicate,  respectively,  that  the  radio  beam  curves 
upward  away  from  the  earth,  downward  with  a curvature  less  than  the 
earth's,  and  downward  with  a curvature  greater  than  the  earth's.  For 
Balboa,  0 < K < 1 occurred  about  2.5  percent  of  the  time,  1 < K < » 
occurred  about  93  percent  of  the  time,  and  K negative  occurred  the  re- 
maining 4.5  percent.  Note  that  the  extreme  K factors  at  Balboa  for 
this  5-year  period  were  +0.39  and  -0.10.  At  Wiesbaden,  0 < K < 1 
occurred  4.6  percent  of  the  time,  1 ^ K < = occurred  95.2  percent  of 
the  time,  and  K was  negative  the  remaining  0.2  percent.  The  extreme  K 
factors  for  Wiesbaden,  +0.36  and  -0.73,  were  not  as  severe  as  those 
found  at  Balboa. 

Equivalent  earth  profile  curves  for  each  of  the  listed  K values  for 
Balboa  and  Wiesbaden  are  also  included  in  figures  4 and  5.  These  sets 
of  profiles  were  constructed  with  the  use  of  equation  11  after  the 
origin  was  chosen,  quite  arbitrarily,  275  ft  above  the  abscissa  which 
was  designated  as  the  flat  earth  surface.  The  coordinates  (d,h)  for 
a selected  K of  interest,  identify  the  location  of  the  radio  terminals. 
Clearance  requirements  above  the  flat  earth  reference  are  usually  dic- 
tated by  the  height  of  obstructing  terrain  and  edifices.  Since  the 
plots  of  figures  4 and  5 represent  the  percentile  deviations  from  the 
median  K values,  these  curves  identify  the  expected  variability  of  radio 
beam  refraction  along  the  path  and  particularly  at  the  radio  terminal. 
The  plots  labeled  +4  and  -4,  accordingly,  represent  the  maximum  earth  - 
beam  divergence  and  maximum  downward  refractivity  which  might  occur. 
Figures  4 and  5 show,  again,  that  radio  beam  dispersion  is  greater  at 
Balboa  than  at  Wiesbaden.  For  example,  if  a transmission  path  of  32 
miles  is  selected,  the  radio  beam  axis  will  remain  within  approximately 
+15  ft  and  -20  ft  of  the  median  K path  63.3  percent  of  the  time  and  +60 
ft  and  -55  ft  95.4  percent  of  the  time  at  Wiesbaden,  Germany,  while  at 
Balboa,  Panama,  the  beam  axis  varies  about  the  median  K path  +20  ft  and 
-35  ft  63.3  percent  of  the  time  and  +70  ft  and  -155  ft,  95.4  percent  of 
the  time.  These  radio  beam  variability  graphs  may  also  be  used  to 
graphically  approximate  the  system  downtime  due  to  atmospheric  refrac- 
tivity once  the  path  length,  antenna  size  and  tower  height,  beam  width, 
and  signal-to-noise  ratio  are  specified. 


EFFECTIVE  EARTH  RADIUS  FACTOR  DISTRIBUTION  FUNCTIONS 

Design  of  high  reliability  microwave  circuits  requires  not  only  a knowl- 
edge of  the  median  effective  earth  radius  factor,  K,  but  also  some  mea- 
sure of  its  distribution  and  the  degree  of  confidence  associated  with 
each  distribution.  Thus,  cumulative  distribution  functions  for  the 
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Figure  4.  Profiles  of  combined  earth  and  refractivity  curvatures  based  on  K median  and  percentile 
deviation  values  derived  from  Balboa,  Panama,  meteorological  data. 
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Figure  5.  Profiles  of  combined  earth  and  refractivity  curvatures  based  on  K median  and  percentile 
deviation  values  derived  from  Wiesbaden,  Germany,  meteorological  data. 


Balboa  and  Wiesbaden  effective  earth  radius  factors  were  constructed 
along  with  the  derived  90  and  99  percent  confidence  interval  for  each 
function.  K factor  cumulative  distributions  were  ordered  from  the 
smallest  positive  K value  to  infinity,  abruptly  changed  to  minus  in- 
finity, and  then  ordered  towards  the  smallest  negative  K,  all  as  a 
function  of  the  percent  of  total  cases.  This  arrangement  is  plotted 
in  the  center  graphs  of  figures  6 and  7,  and  although  it  differs  dras- 
tically from  the  conventional  smallest  to  largest  number  cumulative 
distributions,  it  does  depict  the  degree  and  direction  of  actual  radio 
beam  bending  without  introducing  physical  discontinuities.  Near- 
infinite K values  occur  a significant  portion  of  the  time.  This  ar- 
rangement also  places  the  median  K value  near  the  graph  center  and 
relegates  the  rarer  extreme  refractivity  cases  to  the  plot  edges. 

Figures  6 and  7 are  not  all-inclusive  because  the  few  extreme  valued 
K factors  encountered  called  for  an  unnecessarily  large  graph.  Instead, 
these  K factors  and  the  respective  cumulative  percentages  are  listed  in 
table  2. 

By  selection  of  a specific  Kp  magnitude  from  the  cumulative  distribu- 
tion function,  one  can  predict  the  percentage  of  time  that  radio  beam 
bending  will  remain  within  a required  design  criterion.  However,  each 
Kp  percentile  is  a statistic  based  on  a finite  data  sample.  A small  prob- 
ability exists  that  the  true  distribution  differs  from  the  observed  dis- 
tribution. Thus  confidence  intervals  were  established  to  determine 
within  which  percentile  limits  the  true  Kp  value  would  be  located  given 

a desired  degree  of  confidence.  This  same  Kp  confidence  interval  de- 
fines the  range  within  which  Kp  will  occur  for  the  stated  confidence 

level.  Two  confidence  levels  of  interest  to  USACC  are  90  and  99  percent. 
The  outermost  plots  of  figures  6 and  7 are  the  99  percent  confidence 
limit  boundaries,  while  the  intermediate  plots  are  the  90  percent  con- 
fidence interval  boundaries.  The  interval  widths  are  based  on  the  fol- 
lowing nonparametric  statistical  technique  [7]. 


If  K-j  Kg  . . . Kn  represent  an  ordered  sample  of  independent  random 
variables  with 


K-j  ^ Kg  < . . . ^ Kr  < » . . < Kg  ^ . . . Kn 


and  1 <.  r <_  s ± n;  then  the  lower  confidence  limit  is  given  by 


W a/2  /np(l  - p) 
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TABLE  2.  EXTREME  EFFECTIVE  EARTH  RADIUS  K FACTORS 


Balboa, 

Panama 

Wiesbaden, 

Germany 

Cumulative 

Cumulative 

Percent* 

K 

Percent* 

0.02 

-0.10 

0.02 

-0 

0.04 

-0.17 

0.04 

-0 

0.06 

-0.18 

0.06 

-2 

0.08 

-0.18 

0.08 

-3 

0.10 

-0.25 

0.10 

-3 

0.12 

-0.26 

99.88 

0 

0.14 

-0.33 

99.90 

0 

0.16 

-0.37 

99.92 

0 

0.18 

-0.41 

99.94 

0 

0.20 

-0.42 

99.96 

0 

0.22 

-0.43 

99.98 

0 

0.24 

-0.57 

100.00 

0 

99.87 

-0.58 

99.90 

0.55 

99.92 

0.55 

99.94 

0.54 

99.96 

0.48 

99.98 

0.48 

100.00 

0.39 

* Cumulative  percentages  are  based  on  data  listed  here  and  that  shown  in 
figures  5 and  6. 


Figure  6.  Effective  earth  radius  factor  cumulative  distribution  function  and  confidence 
intervals  for  Balboa,  Panama. 
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and  the  upper  confidence  limit  is  given  by 


s ■ nP  + M1-a/2  /np 


where 


n = total  sample  number, 

p = data  percentile  value  for  which  confidence  interval  is  being 
establ ished, 

a = level  of  significance  or  1 - confidence  coefficient, 

W = the  a/2  or  1 - a/2  percentile  of  a standard  normal  random  variable. 

In  accordance  with  the  above  definitions,  K and  K bound  the  K confidence 

s r p 

interval  for  a stated  1 - a confidence  level. 


Examination  of  the  cumulative  distribution  curves  in  figures  6 and  7 again 
reveals  a higher  variability  in  effective  earth  radius  factor  at  Balboa. 

At  Wiesbaden,  99.86  percent  of  the  effective  earth  radius  factors  encountered 
can  be  accounted  for  by  considering  the  region  0.5  < K <_  °°  plus  -°°  <_  K < -2, 
while  at  Balboa  the  larger  numerical  region  of  0.55  < K <_  °°  plus  -»  <_  K < 

0.58  must  be  considered  to  account  for  99.7  percent  of  the  cases.  The 
Balboa  distribution  function  is  also  more  widely  dispersed  towards  nega- 
tive K values.  Wiesbaden  data  contain  only  a few  negative  K magnitudes, 
which  means  that  there  is  an  infrequent  occurrence  of  severe  downward 
radio  beam  refraction  with  curvatures  greater  than  the  earth's.  Upward 
beam  bending  occurs  approximately  2 percent  of  the  time  at  both  Balboa 
and  Wiesbaden;  this  occurrence  is  demonstrated  by  that  section  of  the 
cumulative  distribution  function  where  0 < K < 1.  Table  2 indicates 
that  the  extreme  downward  refraction  case  for  Wiesbaden,  Germany,  oc- 
curred at  a K value  of  -0.73  which  is  less  severe  than  the  extreme 
Balboa  case  of  K = -0.1.  Table  2 also  shows  that  extreme  upward  beam 
refraction  encountered  at  Wiesbaden  was  for  an  effective  earth  radius 
factor  of  0.36,  while  at  Balboa  the  maximum  upward  refraction  occurred 
at  a comparable  K value  of  0.39.  All  the  extreme  K values  encountered 
represent  refractivity  conditions  vastly  different  from  those  repre- 
sented by  the  respective  medians  or  the  conventionally  used  K value  of 
4/3.  The  variability  of  K is  high,  particularly  at  Balboa,  Panama. 

The  99  percent  confidence  intervals  ranged  from  a AK  of  less  than  0.1 
at  the  median  to  a aK  of  approximately  0.6  at  K = -0.89  (last  interval 


18 


shown  on  figure  6 graph)  for  Balboa,  and  from  a AK  of  approximately  0.1 
at  the  median  to  a AK  of  about  4.3  at  K = -2.0  for  Wiesbaden.  Since  ef- 
fective earth  radius  factor  median  changes  less  than  0.1  are  considered 
negligible  for  present  average  microwave  communications  path  lengths,  a 
99  percent  level  of  confidence  can  be  attached  to  the  fact  that  the  ac- 
tual medians,  both  at  Balboa  and  Wiesbaden,  will  remain  ±0.05  K of  the 
calculated  median.  In  contrast,  the  confidence  intervals  at  the  fringe 
K values  of  -0.89  and  -2.0  can  be  an  appreciable  percentage  of  the  mea- 
sured percentile.  At  Balboa,  the  negative  section  of  the  confidence  in- 
terval at  K = -0.89  is  the  larger  one.  The  magnitude  is  +0.58  and  is 
35  percent  of  the  absolute  value  of  -0.89.  Similarly  at  Wiesbaden  the 
larger  section  of  the  confidence  interval  has  a magnitude  of  approxi- 
mately -3  or  150  percent  of  -2.  The  significance  of  these  percentages 
is  dependent  on  the  circuit  design  criteria  which  define  the  necessity 
to  keep  the  refractivity  within  certain  bounds.  Confidence  intervals 
for  the  extreme  K values  at  Wiesbaden  are  noted  to  be  wider  than  those 
at  Panama.  The  relative  scarcity  or  infrequent  occurrence  of  extreme 
K values  at  Wiesbaden  caused  the  confidence  interval  to  diverge.  Equa- 
tions 12  and  13  indicate  that  the  width  of  the  interval  is  dependent  on 
the  location  of  Kr  and  K$  percentile  points.  The  locations  of  the  Kr  and 

«s  percentiles  are  dependent  on  the  number  of  data  points  available  in 
the  Kp  vicinity.  Consequently,  at  Panama  where  the  higher  K variabil- 
ity places  more  data  points  at  the  fringe  area,  a smaller  interval  de- 
fines the  90  and  99  percent  limits.  At  Wiesbaden  a larger  interval  is 
required  to  retain  the  same  confidence  levels.  Confidence  intervals 
for  the  effective  earth  radius  factor  data  listed  in  table  2 were  not 
derived  because  the  limited  data  sample  in  this  data  group  rendered 
meaningless  statistical  results. 


CONCLUSIONS 


Based  on  an  Army  requirement  to  better  understand  and  quantitatively 
predict  the  degree  of  effective  earth  radius  factor,  K,  variability, 
meteorological  data  were  acquired  from  two  distinct  locations  and  used 
to  derive  the  K factors,  which  served  as  test  cases.  Two  areas  of  con- 
cern to  USACC-EME0,  Germany, and  Panama,  were  chosen  as  the  sample  sites. 
Specifically,  atmospheric  pressure,  temperature,  and  relative  humidity 
data  from  Balboa,  Panama,  and  Wiesbaden,  Germany,  were  converted  to 
average  refractivity  gradients  over  the  first  100  m altitude  with  the 
use  of  an  interpolation  scheme.  Thereafter  for  each  refractivity  gra- 
dient, a corresponding  effective  earth  radius  factor  was  computed.  The 
median  refractivity  gradients  were  determined  to  be  -38  units/km  for 
Wiesbaden  and  -64  units/km  for  Balboa.  The  refractivity  gradient  dis- 
tribution was  found  to  be  non-Gaussian  in  both  cases.  Ninety-five  per- 
cent of  the  Balboa  refractivity  gradients  are  contained  within  -202 
units/km  and  -1  unit/km,  while  at  Wiesbaden  95  percent  of  the  data  is 
within  -86  units/km  and  +16  units/km.  Median  effective  earth  radius 
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factors  were  found  to  be  1.32  for  Wiesbaden  and  1.68  for  Balboa.  Al- 
though the  Wiesbaden  K median  is  in  general  agreement  with  the  current 
"standard"  K = 4/3  magnitude,  the  Balboa  K median  is  significantly  dif- 
ferent. Median  K values  for  84  other  world  locations,  listed  in  table 
1,  demonstrate  the  K factor  median  variability  throughout  the  world.  , 

Derived  K factors  enclosing  68.3,  95.4,  99.7,  and  100  percent  of  the 
Balboa  and  Wiesbaden  K data,  listed  in  figures  4 and  5 along  with  the 
corresponding  equivalent  earth  profiles  curves,  illustrate  the  high 
variability  in  microwave  refraction  at  these  two  locations.  At 
Wiesbaden  radio  beam  upward  refraction  was  found  to  occur  4.6  percent 
of  the  time,  while  downward  beam  bending  occurred  95.4  percent.  At 
Balboa  upward  refraction  occurred  only  2.5  percent  of  the  time,  while 
downward  refraction  occurred  the  other  97.5  percent  of  the  time.  Down- 
ward beam  refraction  is  found  to  be  more  severe  at  Balboa.  K factor 
cumulative  distribution  functions,  shown  in  figures  6 and  7,  also  demon- 
strate the  K variability  at  both  locations.  Boundaries  defining  the 
90  and  99  percent  confidence  limits  of  the  cumulative  distributions 
established  that  both  medians  can  be  predicted  to  occur  within  K < 0.1 
with  a 99  percent  level  of  confidence.  Confidence  intervals  for  extreme 
K values  are  considerably  wider  than  0.1  K.  In  general  K was  found  to 
be  highly  variable  at  both  Wiesbaden  and  Balboa. 

Utilization  of  actual  K medians  at  individual  locations  instead  of  the 
"standard"  4/3  value  is  recommended  when  paths  are  being  established 
for  microwave  circuits.  Such  a change  should  enhance  circuit  reliabil- 
ity. Equivalent  earth  profile  curves  for  various  possible  K values  are 
provided  for  use  as  a test  design  tool.  K factor  cumulative  distribu-  »] 

tion  functions  with  two  confidence  interval  boundaries  are  included  as 
a decision  tool  when  high  reliability  circuits  are  being  prepared. 


* 
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